Evolution of superconductivity in Fe-based systems with doping 
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We study the symmetry and the structure of the gap in Fe-based superconductors by decomposing 
the pairing interaction obtained in the RPA into s- and d-wave components and into contributions 
from scattering between different Fermi surfaces. We show that each interaction is well approximated 
by the lowest angular harmonics and use this simplification to analyze the origin of the attraction 
in and d„^2_y2 channels, the competition between s- and d-wave solutions, and the origin of 
superconductivity in heavily doped systems, when only electron or only hole pockets are present. 
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Introduction. The symmetry and structure of the su- 
perconducting gap in Fe-based superconductors (FeSC), 
and their evolution and possible change with doping are 
currently subjects of intensive debates in the condensed 
matter community. The vast majority of researchers be- 
lieve that superconductivity in FeSC is of electronic ori- 
gin and results from the screened Coulomb interaction, 
enhanced at particular momenta due to strong magnetic 
fluctuations. In FeSC, the Fermi surface (FS) has multi- 
ple sheets due to hybridization of all five d-orbitals of Fe, 
and the interactions between low-energy fermions are a 
complex mixture of contributions from intra - and inter- 
orbital interactions. In this situation, both s-wave and 
non-s-wave pairing are possible, and can be either con- 
ventional or extended, with tt phase shifts between FSs 

Previous theoretical works on FeSCs with hole and 
electron pockets have shown [1" 13| that the s-wave pair- 



ing channel (for sign-changing s gap) is generally the 
most attractive, although the d-wave channel is a strong 
competitor, s-wave gap symmetry is consistent with 
ARPES data, which detected only a small variation of 
the gap along the hole FSs, centered at (0,0), and as 
such ruled out d-wave gap symmetry However, for the re- 
cently discovered heavily electron-doped AFe2Se2(A=K, 
Rb, Cs) fl^ , in which only electron FSs remain according 
to ARPES [il], RPA and functional RG (fRG) studies 
found that the leadi ng p airing instability is now in the 
d-wave channel. 1^ 17 1 d-wavc pairing was also found 
in an fRG study of heavily hole-doped KFC2AS2 [18|, in 
which only hole FSs arc present [l9[. For this material, 
various experimental probes indicate the presence of 
gap nodes, consistent with a d-wave gap symmetry. 

In this communication, we analyze the competition be- 
tween s- and d-wave pairings in doped FeSCs, the origin 
of attraction at small and large dopings, and the struc- 
ture of s- and d-wave gaps at various dopings. We argue 
that the pairing mechanisms at small and large dopings 
are qualitatively different and that the d-wave state at 



large hole doping is a different eigcnstate from the one 
that competes with s-wave at smaller dopings. 

We assume, as in earlier works, that FeSCs can be 
treated as itinerant systems, and that the pairing in- 
teraction is enhanced by spin fluctuations (SF). In the 
band description adopted here, the electronic structure 
at low energies is obtained by hybridization of all five Fe 
d-orbitals and in electron-doped FeSCs consists of two 
cylindrical hole FSs hi and /12, centered at (0,0), and 
two cylindrical electron FSs ei and 62, centered at (7r,0) 
and (0,7r), respectively, in the 1-Fe zone. For hole-doped 
FeSCs, there exists an additional cylindrical hole FS 
centered at (tTjTt). In such a description, interactions 
are dressed by matrix elements associated with the hy- 
bridization of orbitals, and depend on the angles along 
the FSs. 

The method. The input for our analysis is the band 
model with the interactions between the particles on the 
FSs r(ki.,-ki.;k^,-k^) EE r(kF,k^). The interac- 
tions are obtained numerically in the RPA SF formalism 
starting from the 5-orbital model Q with intra- and inter- 
orbital hoppings and local density-density and exchange 
interactions U, U' , J, and J'. We show that, in the band 
basis, each interaction component ry(ki?, k'^;,) is well ap- 
proximated by the leading angular harmonics (LAH) in 
s-wave and d^^-yi-wave channels (similar to the approx- 
imation of the dx2_y2 gap by cos 29 in the cuprates), and 
use the LAH approximation (LAHA) to reduce s-wave 
and d-wave gap equations to either 4 x 4 or 5 x 5 sets 
which can be easily solved and analyzed. This allows us 
to go a step further than previous works and understand 
the pairing mechanism at different dopings, the origin of 
the transition from s-wave to d-wave instability, the role 
of the SF component of the interaction, and the stability 
of s-wave and d-wave gap structures with respect to the 
variation of parameters in the gap equations. For sim- 
plicity, we assume in LAHA that all FSs are circular, with 
the same density of states Np. The results change only 
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a little if we use the actual lattice fermionic dispersion. 

The application of LAHA for FcSCs requires some care, 
as electron FSs are centered at (0, tt) and (tt, 0) points, 
which are not <r^ zLky symmetric. As a result, some 
of the s-wave gap functions, like cos + cos ky behave 
as ± cos 26 along the electron FSs, while some of the d- 
wave gap functions like cos k^ — cos ky are approximated 
on these FSs by constants of opposite sign. With this in 
mind, we treated the angle- independent and cos 2(j) terms 
on equal footings in both s-wave and d-wave components 
of the interactions. A simple analysis then shows that 
LAHA consistent with the FS geometry of FeSCs approx- 
imates the s and dx-2_y2 components of Fi., = NpTij as 
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= U, 



Uhih, + Uh,h, COS 20i cos 2(j)j (1) 

U/i,e(l + 2ah,e COS 26*1) 
+?i/i.e(l + 2afiie cos26'i) cos20i 

(1 -I- 2aee(cOS 26*1 + COS 26*2) + 
4/3ee COS 26*1 COS 26*2) -I- Uee (1 + 
2(5ee(cOS 26*1 -I- COS 2^2) + 4:f3ee COS 29i COS 2^2 



where Uij and Uij are dimensionless interactions in s-wave 
and d-wave channels, respectively, and (f>i and 0i label the 
angles along the hole and electron FSs, measured from 
the k^-axis. Interactions involving other electron FSs are 
obtained by transformations consistent with s-wave or d- 
wave symmetry. 

We use Eq. ^ to fit the RPA interaction r„ by LAHA 
and substitute the parameters extracted from the fit into 
s-wave and d-wavc BCS gap equations, which within 
LAHA are 4x4 matrix equations for two hole and two 
electron FSs and 5x5 when the additional hole FS is 
present. We find the gap structure for the largest positive 
eigenvalue A^.^ (if it exists) and then vary the parameters 
Uij by hand to understand what is the mechanism for the 
attraction. For two hole and two electron FSs the generic 
gap structure is 



A:^(0) = A: + A:cos20, Al_{e) 
At(</')=Af,^cos20, AU^)^A 
Af^(0) = Af + A^os20, Ai^{e) 



Af 



(2) 
Af cos 26* 



h2 COS 20 

: - Af -I- Af COS26I 



and for five FSs we add one more A*J^ {(j)) = A*J^ . 

In Figs. [Hand Owe compare LAHA with the fuU RPA 
Tij{\!iF^'k'p)- The agreement is remarkably good. We 
analyzed eight different sets of U, U' and J, and the 
agreement is equally good for all sets 121[. A very few 
disagreements are cured by adding cos 40 harmonics to 
LAHA. Some of the LAHA parameters extracted from 
the fit, which we will need for comparisons, arc shown in 
Tables U and n 

The cases of weak electron and hole dopings were 
solved numerically within RPA in earlier works, and we 
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FIG. 1; Electron doping, (a-c) Representative LAHA fit of 
tlie interactions r(kF,k^) and a- and d-wave gap functions 
for the case of two very tiny hole pockets. kF is taken to 
be along x on the /12 FS, while k^ is varied along each of 
FSs. The symbols represent the RPA interactions computed 
numerically for the 5-band model using the LDA band 
structure [23|, the black lines are the fits using Eq. [T] Angle 
is measured relative to k^- (d)-(f) are the same as (a)-(c) but 
for stronger electron doping, where there are no hole pockets. 
The parameters are presented in [2TI |. 



TABLE I: Some of the LAHA parameters extracted from the 
fit in Fig. [T] for electron doping. Block (i) corresponds to 
panels (a)-(c) (tiny hole pockets), block (ii) corresponds to 
panels (d)-(f) (no hole pockets). 



(i) 



(ii) 



s-wave Uh-i^hx Uh^e ah^e n 
0.75 0.67 -0.19 0, 



ee <-*ee 



Xs 



Xs 



0.1 0.21 0.84 0.09 -0.12 



d-wave Uh^hi Uhie &hie 



X, 



0.51 -0.32 -0.50 -0.05 0.9 0.35 -0.04 0.1 



Ad 
0.13 



verified that LAHA results are very close to the full 
solutions. For brevity, we present only the results for 
larger dopings, when one type of pockets cither almost 
or completely disappears. We will see that there are quite 
abrupt changes between the two regimes. 

Results and discussion. We varied the magnitudes 
and angle dependencies of the interactions by hand and 
checked what most influences the value of A and the 
structure of the gap. We found that some system prop- 
erties are sensitive to the ratios of the parameters, but 
some are quite universal. 
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For electron doping, parameter-sensitive properties in- 
clude the gap symmetry, since As and remain compa- 
rable as long as both hole and electron FSs are present 
(see Table m, and the presence or absence of accidental 
nodes in the s-wave gap, although for most of parameters 
the gap does have nodes, as in Fig. [TJb). The universal 
observation is that the driving force for attraction in both 
s-wave and d-wave channels is the inter-pocket electron- 
hole interaction (u^.g and u^.g terms), no matter how 
small the hole pockets are. When the SF component of 
the interaction is large, Uh^e and Uhic exceed the hole- 
hole and electron-electron interactions. Then X^^d are 
positive already if we neglect the cos 26 terms in ^ (for 
two equal hole FSs the conditions are > uthUee and 
> u/ihUee)- In this case, the cos 20 terms in the s- 
wave and d-wave gaps scale with the corresponding ahe- 
For smaller SF component, when uf^^ < UhhUee (the case 
considered in Fig. [T] and Table H]), the electron- hole inter- 
action still generates solutions with Xg^d > 0, only this 
time the gap develops a stronger cos 29 component, which 
effectively reduces Uee- 

The situation changes qualitatively once the hole pock- 
ets disappear (Fig. [T]Jd)-(f)). We see from Table |T] that 
Xs is reduced, but Ad is enhanced, i.e., the c?-wave Tc 
increases. Comparing the LAHA parameters for the 
two dopings, we see the reason: once the hole pock- 
ets disappear, a direct d-wave electron-electron interac- 
tion Uee becomes strong and attractive. To understand 
why this happens, we note that Uee and Uee are sym- 
metric and antisymmetric combinations of intra-pocket 
and inter-pocket electron-electron interactions: Uge = 



Both 



and 



are positive (repulsive), hence Uee > 0, but the sign of 
Uee depends on the interplay between uf^^.^J. and uf^^.J.^. 
As long as the hole FS is present, SF are peaked near 
q = (0,7r) and (7r,0), which are an equal distance from 
the relevant momenta q = for uf^^^.^ and q = (tt, tt) for 
""fntor- ^his situation, u'?^^^^ and u?^^.^J. remain close in 
magnitude, and Uee is small. Once the hole pocket dis- 
appears, the peak in the RPA spin susceptibility shifts 
towards (tt, tt) and uf^^.^J. increases more due to the 
SF component than uf^tra- ^ negative uf^^.^^ — uf^j^r then 
gives rise to a "plus-minus" gap on the two electron FSs 
for the same reason that large Uhe gives rise to a sign- 
changing gap between hole and electron FSs, (and the 
interaction between hot spots in the cuprates gives rise 
to a sign-changing gap in the hot regions). Such a gap 
changes sign between electron pockets, which diffe r by 
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kx ky and therefore has d^^-yi symmetry [1 
Our gap functions in both s- and d-wave channels at large 
electron doping are in good quantitative agreement with 
the full solution of the RPA gap equation [l3| and with 
fRG results [H. 

Next we consider the case of hole doping. The LAHA 
fits to the cases when electron FSs are small but still 
present and when only hole FSs remain are shown in 
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FIG. 2: The same as in Fig. [U but for hole doping (3 hole 
FSs). Panels (a)-(c) are for the case of tiny electron pockets, 
(d)-(f) are for stronger hole doping, when there are no electron 
pockets. The parameters are presented in [2ll ]. 



TABLE II: Some of LAHA parameters extracted from the fit 
in Fig. [2] for hole doping. Block (i) corresponds to panels (a)- 

(c) (tiny electron pockets), block (ii) corresponds to panels 

(d) -(f) (no electron pockets). 



(i) 



S Uhihi Uhie Cthie Uee \a 

0.75 1.36 0.08 1.40 1.8 

d Uhihi Uhie &hie Uee Ad 

0.70 -1.32 0.0 1.45 1.2 



(ii) 



Uhihi Uhih.2 Uhiha Uh^h3 Xs 

0.67 0.8 0.29 1.37 0.13 

Uhihi Uhih2 Uhih3 Uh-ih^ Xd 

0.36 -0.5 -0.02 -0.17 0.11 



Fig. [2l The parameters extracted from the fit are shown 
in Table |IT1 We analyzed these and other dopings and 
again found universal and parameter-sensitive features. 
The parameter-sensitive property is again the presence 
or absence of accidental nodes in the s-wave gap along 
the electron FSs. For most of the parameters, the gap 
does not have nodes (see Fig. [5]) because Uhe increases 
once it acquires an additional contribution it/ige, but for 
some parameters we still found nodes along the electron 
FSs. The universal observations are that, as long as both 
hole and electron pockets are present, (i) the s-wave is 
the leading instability (As > A^ > 0), and (ii) the driv- 
ing force for the attraction in both s- and d- channels is 
again inter-pocket electron-hole interaction (uhe and Uhe 
terms), no matter how small the electron pockets are. In 
the d-wave channel, the electron-hole interaction changes 
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sign between the two hole FSs at (0, 0), as a result d-wave 
gaps on these FS have a 7r-phase shift (see Fig. [2l^c)). 

The situation rapidly changes once electron pockets 
disappear. The d-wave eigenvalue grows relative to 
As and for the doping shown in Fig. [5] almost exceeds it. 
It is very likely that d-wave becomes the leading instabil- 
ity at even higher dopings, and we therefore focus on the 
d-wave channel. Comparing u in Table HIl for the cases 
with and without electron pockets, we find that the d- 
wave channel is attractive in the absence of the electron- 
hole interaction because of two reasons. First, the d-wave 
intra-pocket interaction Uh^h^ becomes negative (attrac- 
tive). Second, the inter-pocket interaction u^i/ia larger 
in magnitude than repulsive Uh^hi and Uh2h2- The solu- 
tions with positive A^ then exist separately for FSs hi^2 
and /13, and the residual inter-pocket interaction just sets 
the relative magnitudes and phases between the gaps at 
/13 and /ii_2. Because Uhi^h2 is attractive, the two d-wave 
gaps at are now in phase, i.e., this d-wave solution 
is a different eigenfunction from the one with phase shift 
TT at smaller dopings. The difference is clearly seen by 
comparing panels (c) and (f ) in Fig. [31 The d-wave gap 
symmetry at large doping and in-phase structure of the 
gaps at hi ^2 is consistent with the fRG solution [3 

Conclusions. The key result of this work is the obser- 
vation that the mechanism of the pairing in FeSCs with 
hole and electron FSs is different from the one at strong 
hole or electron doping, when only one type of FS re- 
mains. At small/moderate dopings, the pairing is driven 
by inter-pocket electron-hole interaction, no matter how 
small hole or electron FSs are. In hole-doped FeSCs, 
the leading instability is s-wavc, while in electron-doped 
FeSCs, s- and d-wave channels are strong competitors, 
and which of the two wins depends on the model pa- 
rameters. At large electron and hole dopings, d-wave is 
the leading instability, although the s-wave channel re- 
mains attractive. At strong electron doping, the origin 
of the pairing is a direct d-wave attraction between elec- 
tron pockets. At strong hole doping, however, the rea- 
son for the d-wave pairing is a d-wave attraction within 
the (tt, tt) pocket and between the two hole pockets at 
(0,0). The d-wave pairing at strong hole doping is con- 
sistent with the observation of nodal quasiparticles [2^ 
in the heavily hole doped superconductor KFe2As2with 
Tc = 3K. Superconductivity at heavy electron doping at 
a rather high Tc 30K has been recently discovered in 
AFc2Se2(A:=K, Cs. Rb), which onl y h ave electron FSs, 
according to recent ARPES studies [l5| . Whether this is 
a d-wave superconductor remains to be seen. 

We have only studied the strictly 2D case thus far, and 
neglected aspects of the 3D I4/mmm crystal symmetry 
characteristic of 122 materials and the hybridization of 
electron pockets in the folded zone [i^] . We nevertheless 
believe that the general evolution of interactions and gap 
symmetry discussed here will be generic to the FeSCs. 
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